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ELASTIC SCATTERING OF SLOW ELECTRONS 
BY TWO-ELECTRON IONS 

M. R. C. McDowell 
ABSTRACT 

A one channel model, previously used for He, is extended to allow calcula- 
tion of elastic scattering by H and Li^ at energies up to three rydbergs. Non- 
separable ground state wave functions are used, allowance beii^ made for ex- 
change, and dipole and quadrupole polarization potentials. Results for Li^ are 
in good agreement with quantum defect method values. Calculated differential 
cross sections for elastic scattering by H at energies of one half and one 
rydberg are presented. Deviations from Couloumb scattering are marked. 
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ELASTIC SCATTERING OF SLOW ELECTRONS 


BY TWO-ELECTRON IONS 

by 

M. R. C. McDowell*^ 
NASA-Goddard Space Flight Center 
Greenbelt, Maryland 


I. Introduction 

A considerable number of theoretical investigations on elastic scattering 
of slow electrons by He atoms have been reported recently. They are in 
substantial agreement with each other, and are consistent (in the sense of dis- 
persion relations^) with experimental data.5-6 In view of advances in experi- 
mental technique which have allowed measurement ^ ® of inelastic collision cross 
sections of slow electrons with H (in good agreement with theory it is of in- 
terest to examine elastic scattering by two electron ions. 

In this paper the model of Williamson and McDowell^ has been extended and 
applied to elastic scattering by H and by Li^ at energies up to 3 Ry. The theory 
is outlined in Section II, and the numerical methods employed are discussed in 
Section III. Phase shifts and scattered intensities are presented in Section IV, 


*National Academy of Sciences-National Research Council, Senior Post-doctoral Associate 
^On leave of absence from University of Durham, England 1967-8. 
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and in the case of Li^ compared with the result of a Hartree-Fock calculation^®, 
and the quantum defect method. 

n. Theory 

We consider an electron scattered by a two electron system of charge z. 
The total three electron system is described by the Schrodinger equation, 

(H - E) ¥ (1, 2, 3) = 0 ( 1 ) 

with Hamiltonian operator* 
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In our model the total wave function (123) is represented by an ansatz, 

'P (1,2,3)= '/'o(l-2) F(3) S(123) (3) 

1 , 2,3 

where (1, 2) is a wave function for the ground state of the target when elec- 
trons 1 and 2 are bound, S (123) is a spin function, and F (3) the unknown scat- 
tering function to be determined, and the sum is over cyclic permutations. 


*AII quantities will be given in atomic units, the unit of energy being the Rydberg, 13.565 eV. 
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Taking 


S(123) -a2/3i)a3 (4) 

where a.,/3 are the one-electron spin up and spin down functions, we may write, 
in a usual notation 

F(3) = 

^=0 3 

the position vector of electron 3 with respect to the target nucleus being 
Tg = (r 3 , ^3 , 4>z) and ^3 = cos^g. 


The integrodifferential equation for the scattering functions (r 3 > is ob- 
tained by projecting out ( 1 ) on each partial wave in turn, 

j*0*o ( 1 , 2 ) (/L 63 ) S* (123) [H-E] W(123) dr^ dr^dr^ d Spin = 0. ( 6 ) 

For simplicity we assume that the ground state function for the target is known 

I 

exactly: 


(- 




^2 12/ 


and 


E = E. +k2. 


We adopt a two-parameter variational trial function 

,/;,(1.2) =N (e 


-arj -/3f2 ^ 


+ e 


(7) 


( 8 ) 


(9) 
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for the target. Then on carrying out the angular and spin integrations (6) yields 




16 77 ^ N ^ 

(2't + l) 



fo(l)Zoo(123)dr, dr^ 


( 10 ) 


n co 

f^(l)Xo^(123)d. 

) 


dr. 




where 


is the operator 


+k2 

dr2 r2 


-Voo (3)- 


( 11 ) 


The direct and exchange potentials are 


Voo (O = - 


2z 


4t7^ 

a^/S^T 


1 - (1 + ar] + [1 - (1 + /3 r) 


^1280^ [i _ fi +ir;7r)e-<“+^)' 
(a+/3)6 L \ 2 J 


^00 “ 


Ta ^ ^ -cir.-jSr^ ^ -/3rj-ar 

^ (23) r2 J(a2rj - 2 a) e ^ r ^ - 2/5) e 


^*<,(1,2) j^r,kt +2z-^ j-, 


X„t = t,r,r=0;(l,2).jA- 'A„(2.3), r,,^; 



( 12 ) 
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No account has been taken of the polarizability of the target in this model, 
so we replace ^ by 

■ 2 (Vd + \) (13) 

where Vd and are potentials behaving as r"'* and r~® respectively at large 
r. They are chosen in the Bethe-Reeh form^^ 



Vq (r) = 


-iL(l-2ye-v r 
2y6 I L 


2y e 2y I 1 + 2y + 


+ — y® + — y® - y^ + — ^ y® - y® 
9 ^ 18 175 135 135 


] 


_4vP« M 20 2 c 3 284 85 Ig"! 

e M + 4y +— y^ + 6y^ ^ '^ 9 '^ '9^1 
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where Ei (-u) = -E^ (u)and Ei(u) is the exponential integral, 

u = Zj r , y = Zj r (15) 

The parameters z^, are chosen so that (9/zj)^ =°‘d> (IS/zj)® = a^, where 
a^, are the dipole and quadrupole polarizabilities of the target respectively. 

III. Numerical Methods 

The integrodifferential equation 

must be solved subject to the boundary conditions 

(0) = 0, (r) ~ sin ^kr + "j“ 2kr ~ ^ 

with 

<^j^ = avgl^ (z + 1--^ 

and Zq = (z - 2). The quantity of interest is the non-couloumb part of the phase 
shift, 5^ . We shall refer to it simply as the phase shift (for the partial 
wave). 

The equation (10') was solved by a non-iterative method. Writing 

(r) = P (r) + /i-Q (r) + vR (r) (17) 
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where P, Q, R satisfy 

L^P = F(P) 

L^Q =F(Q) + 

L^R=F(E)tc^r^*'e-/>',c^ = |^ ( 18 ) 

and F (0) involves linear combinations of integrals over the range (0,r) only, 
the parameters ft, v may be determined in terms of certain infinite integrals 
once P(r), Q(r) and R(r) are known. For He this procedure yields phase shifts 
identical to those obtained previously^ by an iterative method. The ordinary 
differential equations (18) were solved by a Fox-Goodwin predictor-corrector 
method. 

The solutions were normalized by the Stromgren procedure^'* as given by 
Burgess. His procedure for determining the phase shift when a polarization 
potential is present may be extended easily to the case of a negative ion (zg = -1). 
At large r, 

k"^/^sin[c/)(r)+Sp] (19) 

^ r-»oo ^ 

where to sufficient accuracy 

0(r) = (r) +1 <E>(r). 
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Defining 


yo = lz(,r I , /3 = |z(jkl 
we have (the ± sign being that of z^), 

=X(^)* - +^X* ± 


with 


^ 5(yO"c) _ X* (3/3^ +4) + /3p (3y5^c +2) ±ygc 
24^(±)3 24x (1 c) (X* +/3p) 



( 20 ) 


X* = + (/3^ ± 2p - , T* (/o) = X. + c +/0 ^ c = -^' ("t + 1) . 

(1 +/?2c) 


The polarization term is 


<D(yO) = 




4c^ 


3(2p;c) ^cx^ 


(3 +/3^c) 


COS 


-1 


,p(x*+/3p) P' 


v^C 


T*(P) 


, c 0 


( 21 ) 


2z^ ttj {9!ipx^ + ll/32p^ ± 6p) 

15 p2 (X* +/3p)3 


, c = 0. 


When c = 0 the cos"^ term in (20) is replaced by 0.25 (x* +y5p)”^- No special 
treatment of the quadrupole potential is required. Dipole and quadrupole polar- 
izabilities for H are given by Stewart.^ The value of the quadrupole polariz- 
ability used is uncertain by ±20%, but in view of its small effect, this is not 
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significant. The dipole polarizability used for Li^ is the Hartree-Fock value of 
Lahiri and Mukkerji^® and should be accurate; its quadrupole polarizability was 
taken as zero. The calculated phase shifts for H are shown in Table 2. The 
program actually calculates tan S^, and (mod 77). The calculations indicated 
that if the phase shift was defined to go to zero in the high energy limit, then the 
s-wave phase shift must be chosen to be 77 in the zero-energy limit. This was 
confirmed by comparing the numerical solution for k^ = 0.05, ■I - 0 with the cor- 
responding pure couloumb solution. An extra node is present in the calculated 
solution at small r . Physically this occurs because the incoming electron 
cannot enter the filled Is-shell, and there is in our model, no bound (Is) 2s 
state of H .A similar situation occurs for electron scattering by He, but not 
for Li . In the latter case the (Is) 2s state is the ground configuration of Li, 
and the zero energy s-wave phase shift is found to be 1.261 (see below). The 
s-wave phase shifts for all these systems are shown in Fig. 1. In Fig. 2 we 
show tan for H both including (d + q) and excluding (d) the quadrupole poten- 
tial. The local potentials occurring in Eq. (10’) are shown in Fig. 3. They are 
strongly repulsive at large r but the effect of the attractive polarization poten- 
tials is to produce a broad positive maximum in V(r) from r = 2.0 to r = 10.0 
with a maximum value of 0.2 Ry. The potential V(r) then becomes strongly 
attractive for r < 2.0 a^^. There appears to be a possibility of trapping for low 
energy electrons. We therefore computed tan for 0.05 < k^ < 0.2 Ry at 
intervals of 0.01 in k^ (i.e. 0.136 eV). No evidence for an s-wave resonance 
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was observed, but in any event no resonance appreciably narrower than 0.1 eV 
would have been found. Herzenberg and Lau^^ have used our model (with a dif- 
ferent form of polarization potential) to investigate s-wave scattering of slow 
electrons in He, and obtain a resonance by suitably modifying the strength of 
the exchange terms. Peterkop^® has also reported an s-wave resonance in He, 
in a one channel model. 

Although the calculated results for He and Li^ are in excellent agreement 
with other calculations (and for He, with experiment), it might appear that the 
model would be less satisfactory for H , with its extremely large polarizability. 
In fact the calculated phase shifts in He and Li are also dominated by the polar- 
izability potential (for small ^ and k^) and are substantially in error when 
polarization is neglected. We would therefore expect our H results to be of 
comparable accuracy. 

The effect of introducing a quadrupole polarization is to slightly increase 
all the phase shifts, in general by less than 10%. The highest value of k^ at 
which the p-wave phase shift passes through tt/ 2 increases from 2.0 to 2.45 
when changes from zero to 1300. For given k^, the phase shifts decrease 
slowly with - 1 ). 

Calculated phase shifts for Li^ are displayed in Table 3. When the polar- 
ization potentials are omitted the calculated values are in very close agreement 
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with these obtained in a Hartree-Fock calculation.^® This confirms the result 


obtained for He^ that provided the target ground state wave function used gives 
a binding energy at least as good as the Hartree-Fock value, its detailed behavior 
is not important in determining the low energy phase shifts. In particular it ap- 
pears that short range -1 = 0 correlation in the target is not important. The Li^ 
p-wave phase shift obtained when polarization is included tends to a threshold 
value (0.170) in close agreement with that (0.188) obtained using the quantum 
defect method, but rather larger than the no-polarization value of 0.110. 

The zero energy s-wave phase shift is 1.261 where the quantum defect 
method gives 1.232 ± 0.005. The d-wave phase shift at threshold, is like its 
quantum defect value, very small (approximately 3.8 x lO”^). 

We have computed angular distributions of elastically scattered electrons 
for H , where deviations from pure Couloumb scattering would be expected to 
be significant. 

The non-couloumb contribution f^(6') to the scattering amplitude f(6>) may 
be written 

00 

0 

and because of the large polarizability of H a comparatively large number of 
terms must be retained in the summation over Z in (22). We find that for 
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- 3.0, the sum has converged to better than 1% if it is cut off at -i = 30. For 
'1^4 the exchange terms in (10') may be ignored; phase shifts are then obtained 
from the suitably normalized solutions of 

L^(r) f^(r) = 0. (23) 

In table (2b) we give phase shifts in this approximation for -t = 4 (1) 10 at 
= 0.5, 0.75, 1.0, 3.0. 

Rather than displaying the scattering intensity 1{6) = lf(6>) 1^ it is convenient 
first to show R((9), the ratio of 1(6') to the pure coloumb intensity (6). Writing 
a = l/k we have^° 

R (6,k2) = 1 1 + n|^ (24) 


with (for H ), 



The calculated values of R(6, k^) for k^ = 0.50, 0.75, 1.0, 3.0 are shown 
in Fig. 4. 


The most significant features are the sharp dip near 6 = 0.4, the oscilla- 
tions in tiie range 1.0 < 6 5 2.5 and the strong enhancement in the backward 
direction for = 0.5. The fine detail of the calculated oscillations in R(6, k^) 
at intermediate angles may not be significant. 
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In Figs. 5a, 5b we show 1{6) and 1^(0) for = 0.5 and 1.0 over the angular 
range 0.3 S 0 ^ tt. The calculated intensity is a factor of five lower than the 
couloumb intensity at 40°, but for = 0.5 it is a factor of four higher in the 
backward direction. 
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Figure Captions 


Fig. 1 Calculated s-wave phase shift for elastic scattering of slow electrons 
by the two electron systems indicated. 

Fig. 2 Tan 8^ (k for elastic scattering by H . Full curve = 0, dashed 

curve a ^ 0 . 
q 

Fig. 3 The local potentials occurring in Eq. (10’) for Z = 1 (H~). 

Fig. 4 Calculated values of R(i9 , k*), the ratio of 1(0) to the Couloumb 

intensity (^) for h", (1) k* = 0.5 (2) k^ = 0.75 (3) k^ = l.O 
(4) k2 = 3.0. 

Fig. 5 Calculated differential cross sections for scattering by H . The full 
curve is 1^(0), the dashed curve 1(0). (a) k^ = 0.5, (b) k^ = 1.0. 



Table 1 


Parameters for H , li^ 


a 

/3 


a 

a 

Z 

Zo 

1.04 

0.2808 

203 

1300 

1 

-1 

3.295 

2.079 

0.19 

0 

3 

+1 
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Table 2 a 


S-Wave Phase Shifts for H 



tan 


^0 

k2 

d 

d+q 

d+q 

0.05 

8.19, -3 

8. 22, -3 

3.150 

0.06 

1.26, -2 

1.27, -2 

3.154 

0.07 

1.84,-2 

1.85, -2 

3.160 

0.08 

2.51, -2 

2.54, -2 

3.167 

0.09 

3.26, -2 

3.31, -2 

3.175 

0.10 

4.06, -2 

4.13, -2 

3.183 

0.25 

2.22, -1 

2.50, -1 

3.389 

0.50 

4.71, -1 

5.59,-l 

3.650 

0.75 

3.49,-l 

4.20,-l 

3.540 

1.0 

1.74, -1 

2.30,-l 

3.369 

2.0 

-3.49, -1 

-3.07,-! 

2.844 

3.0 

-7.43, -1 

-6.97, -1 

2.510 


(d indicates dipole potential only 
d + q indicates dipole plus quadrupole potentials) 
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Table 2b 


Calculated H Phase Shifts 


= 1 


k2 

d 

d +q 

0.10 

4.32, -2 

4.38, -2 

0.55 

1.49, -1 

1.53.-1 

0.20 

3. 72,-1 

4.00,-! 

0.25 

7.28, -1 

8.12,-! 

0.50 

1.71 

1.82 

0.75 

1.78 

1.86 

1.0 

1.75 

1.82 

1.5 

1.65 

1.70 

2.0 

1.57 

1.61 

3.0 

1.44 

1.46 


Notation as in Table 2 a 

^ = 2 


d 

d +q 

2.95, -2 

2. 99, -2 

7. 41, -2 

8. 08, -2 

1.59, -1 

1. 63,-1 

2.64, -1 

2.76, -1 

8.13,-! 

8.78, -1 

1.04 

1.11 

1.11 

1.18 

1.12 

1.17 

1.09 

1.14 

1.03 

1.07 
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Table 2 c 


Dipole Plus Quadrupole Potentials Only 


k2 

= 3 

^ =4 

^ = 4 

no exchange 

0.25 

0.139 

0.0798 

- 

0.5 

0.400 

0.221 

- 

0.75 

0.667 

0.403 

0.390 

1.0 

0.768 

0.503 

0.497 

1.5 

0.858 

0.624 

- 

2.0 

0.869 

0.669 

- 

3.0 

0.846 

0.687 

_ 
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Table 2d 



0.50 

0.75 

1.0 

3.0 


Phase Shifts for H , 'I t 4, no exchange 
(dipole plus quadrupole potential only) 


4 

5 

6 

7 

8 

9 

10 

0.230 

0.147 

0.0944 

0.0637 

0.0449 

0.0327 

0.0244 

0.390 

0.244 

0.158 

0.107 

0.075 

0.0540 

0.0403 

0.497 

0.329 

0.221 

0.152 

0.107 

0.0775 

0.0577 

0.687 

0.552 

0.447 

0.360 

0.288 

0.231 

0.185 
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Table 3 


Phase Shifts for Dipole Potential Only 


k2 

o 

II 

1 = 1 

1 =1,H.F. 

II 

P 

a 

II 

O 

1 = 2 

0 

(1.261) 

(0.170) 

(0.110) 

- 

(0.0038) 

0.04 

1.258 

0.173 

0.113 

- 

0.0044 

0.09 

1.255 

0.177 

0.117 

- 

0.0052 

0.16 

1.248 

0.182 

0.120 

0.122 

0.0063 

0.25 

1.241 

0.189 

0.126 

- 

0.0074 

0.36 

1.230 

0.195 

0.132 

0.133 

0.0088 

0.49 

1.221 

0.202 

0.140 

0.140 

0.0108 

0.64 

1.209 

0.210 

0.147 

0.146 

0.0135 

0.81 

1.194 

0.217 

0.153 

- 

0.0157 

1.0 

1.178 

0.224 

0.161 


0.0182 


(1) Values for = 0 are obtained by graphical extrapolation. 

(2) For the p-wave phase shift the successive columns show 
the results (i) in our model with = 0.19, (ii) in a Hartree- 
Fock calculation,^® (iii) in our model with = 0. 
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Figure 3. The local potentials occurring in Eq. (10') for 












